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ABSTRACT 


The  finite  segment  modelling  of  a  long  flexible  cable  is 
discussed.  An  equivalent  force  system  representing  the  fluid  and 
gravitational  forces  exerted  on  a  typical  cable  segment  is  determined. 
The  fluid  forces  include  inertia,  drag,  and  hydrostatic  (bouyancy) 
forces.  Application  of  the  results  in  computer  algorithm  development 
is  also  discussed. 


\r-.J 


INTRODUCTION 


The  objective  of  this  report  is  to  develop  expressions  for  the 
fluid  forces  on  a  submerged  segment  of  a  finite  segment  cable  model. 

It  is  expected  that  these  expressions,  when  used  with  the  computer 
algorithms  of  a  chain  link  model,  will  provide  a  comprehensive  yet 
efficient  simulation  of  submerged  cable  dynamics. 

There  has  been  an  increasing  interest  recently  in  the  dynamics  of 
submerged  cables.  This  interest  is  due  to  the  significant,  and  often 
deleterious,  effects  of  the  fluid  forces  on  cable  motion  and  cable 
dynamics  in  a  wide  range  of  towing  and  mooring  configurations . 

This  interest,  in  turn,  has  stimulated  the  development  of  a  number  of 
procedures  for  modelling  and  studying  cable  dynamics  for  a  variety 
of  applications.  References  [1-353*  provide  a  brief  summary  of 
some  of  these  procedures.  In  a  review  of  these  procedures,  Choo  and 
Casarella  I10J  suggest  that  a  promising  approach  to  obtaining  an 
accurate  analysis  of  large  displacement,  unsteady,  three-dimensional 
motion,  is  the  finite  segment  (that  is,  rigid  link)  modelling  of  the 
cable. 

Independent  of  these  modelling  procedures  in  cable  dynamics, 
there  have  also  been  many  recent  significant  advances  in  finite  seg¬ 
ment  modelling  of  other  mechanical  systems  --  particularly  "multibody 
systems."  A  multibody  system  or  "chain  system"  (sometimes  also  called 


♦Numbers  in  brackets  refer  to  References  at  the  end  of  the  report. 
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"open-chain"  or  "cpen-tree"  systems)  is  defined  as  a  set  of  rigid 
bodies  arbitrarily  assembled  such  that  adjacent  bodies  share  a  common 
point  and  such  that  no  closed  loops  are  formed.  Figure  1.  presents  a 
sketch  of  such  a  system. 

The  interest  in  general  chain  systems  stems  from  the  fact  that 
they  are  excellent  models  of  many  physical  systems  such  as  manipulators, 
cranes,  robots,  biodynamic  systems  (for  example,  human  body  models), 
chains,  and  cables. 

The  advances  in  the  modelling  and  analysis  of  such  systems  are 
due  to  a  large  extent  to  corresponding  advances  in  computational 
technology  and  in  computer  hardware.  But,  beyond  this  there  have  also 
been  significant  advances  in  the  procedures  for  developing  the 
governing  dynamical  equations  of  motion  of  these  systems.  One  of  the 
most  promising  of  these  procedures  is  based  on  using  Lagrange's  form 
of  d'Alembert's  principle  to  obtain  the  equations  of  motion  (See 
References  136-4(2.).  This  principle,  as  exposited  by  Kane  and  others 
[42-44]  is  similar  to  Lagrange's  equations  in  that  non-working  internal 
constraint  forces  between  the  bodies  of  the  system,  are  automatically 
eliminated  in  the  analysis.  Also,  the  principle  leads  directly  to 
the  correct  number  of  governing  equations  (that  is,  one  equation  for 
each  degree  of  freedom).  However,  unlike  Lagrange's  equations,  the 
principle  does  not  require  the  differentiation  of  lengthy  scalar 
energy  functions.  Instead,  it  uses  vector  derivatives  which  can  be 
calculated  by  vector  multiplication.  Such  calculations  are  ideally 


suited  for  conversion  into  computer  algorithms.  Indeed,  in  References 
[36-40] ,expl icit  expressions  and  algorithms  for  obtaining  the  coefficients 
of  the  governing  dynamical  equations  are  presented. 

This  approach  is  directly  applicable  in  the  modelling  and  analysis 
of  cable  dynamics.  However,  to  develop  a  cable  model,  that  is,  a  link 
or  finite-segment  cable  model,  it  is  necessary  to  have  expressions 
for  the  fluid  forces  acting  on  the  cable  segments.  Hence,  the  emphasis 
of  this  report  is  the  analysis  of  these  fluid  forces  and  the  develop¬ 
ment  of  analytical  expressions  for  them  --  in  a  format  suitable  for 
their  direct  incorporation  into  multi  body  computer  models  as  described 
in  References  [38-40,45,46]. 

The  balance  of  the  report  is  divided  into  six  sections.  The 
first  of  these  provides  the  assumptions  and  the  basis  of  the  analysis 
used  in  the  sequel.  The  next  section  describes  the  fluid  inertia 
forces  transmitted  to  the  cable  and  generated  by  acceleration  of  the 
cable  segments.  This  is  followed  by  two  sections  describing  the  normal 
and  tangential  drag  forces  and  the  weight  and  bouyancy  forces  on  the 
segments.  The  final  section  contains  some  concluding  remarks. 


PRELIMINARY  CONSIDERATIONS 


NOTATION 

a,b,c  -  Coefficients  defined  by  Equations  (32),  (33),  and  (34). 

p 

a  -  Acceleration  of  the  mass  center  G  of  a  cable  segment  in 
an  inertial  reference  frame  R. 

-GN  "  Acceleration  of  the  fluid  relative  to  G  in  R  and  normal  to 
the  cable  segment  axis. 

-  Normal  component  of  the  acceleration  of  the  fluid  relative 
to  the  cable  segment. 

p 

a  -  Acceleration  of  a  typical  point  P  of  a  cable  segment  in  an 
inertial  reference  frame  R. 

W 

a  -  Acceleration  of  the  fluid  in  an  inertial  reference  frame  R. 

aW^G  -  Acceleration  of  the  fluid  relative  to  the  cable  segment  at  G, 

the  mass  center  of  the  segment. 

W/P 

a  -  Acceleration  of  the  fluid  relative  to  the  cable  segment  at 
a  typical  point  P. 

A,B,C  -  Fluid  force  coefficients  of  Equations  (1),  (2),  (3),  and  (4). 

B  -  Bouyancy  force  on  a  submerged  right  circular  cylinder. 

B^  '  -  Component  of  B  normal  to  the  cable  segment  axis. 

By  -  Component  of  B  parallel  to  the  cable  segment  axis. 

CM  -  Added  mass  coefficient. 

CN  -  Normal  drag  coefficient. 

Cy  -  Tangential  drag  coefficient. 

d  -  Diameter  of  the  cylindrical  cable  segment. 

f  -  Fluid  force  per  unit  length  exerted  on  a  cable  segment. 

Fw  -  Resultant  added  mass  force  passing  through  the  mass  center  of 

a  cable  segment. 


F  -  Resultant  fluid  and  gravitational  force  passing  through  the 
mass  center  G  of  the  cable  segment. 

F*.  -  Resultant  normal  drag  force  passing  through  the  mass  center  G 

of  the  cable  segment. 

Ft  -  Resultant  tangential  drag  force  passing  through  the  mass 
center  G  of  the  cable  segment. 

g  -  Gravity  constant. 

G  -  Mass  center  of  a  typical  cable  segment. 

k  -  Vertical  unit  vector. 

L  -  Length  of  a  typical  cable  segment. 

n  -  A  unit  vector  parallel  to  the  cable  segment  axis. 

P  -  A  typical  point  on  the  cable  segment  (See  Figure  3.). 

R  -  An  inertial  reference  frame. 

Rg  -  Normal  Reynolds  number  (See  Equation  (8)). 

ReT  "  Tan9entia1  Reynolds  number  (See  Equation  (9)). 

T  -  Torque  of  the  couple  of  the  equivalent  fluid  and  gravitational 
force  system. 

JM  -  Torque  of  the  couple  of  the  equivalent  added  mass  force  system. 

-  Torque  of  the  couple  of  the  equivalent  normal  drag  force  system. 

T,  -  Torque  of  the  couple  of  the  equivalent  tangential  drag  force 
system. 

V  -  Cable  segment  volume. 

V_  -  Velocity  of  the  mass  center  G  of  a  cable  segment  in  an  inertial 
reference  frame  R. 

VGN  "  Velocity  of  the  fluid  relative  to  G  in  R  and  normal  to  the 
cable  segment  axis. 

VN  -  Normal  component  of  the  velocity  of  the  fluid  relative  to  the 

~n  cable  segment. 

p 

V  -  Velocity  of  a  typical  point  P  of  the  cable  segment  in  an 

inertial  reference  frame  R. 

9*  T  «***•*'  '  ■*  -  - _ 


i  x: 


VP/G  -  Velocity  of  the  fluid  relative  to  the  cable  segment  at  G, 
the  mass  center  of  the  segment. 

V,  -  Tangential  component  of  the  velocity  of  the  fluid  relative 
to  the  cable, 
u 

V  -  Velocity  of  the  fluid  in  an  inertial  reference  frame  R. 

Resultant  weight  force  on  the  cable  segment, 
x  -  Coordinate  of  P  relative  to  G. 

X  -  Defined  by  Equation  (40). 

a  -  Angular  acceleration  of  a  typical  cable  segment, 

u  -  Viscosity  of  the  fluid, 

p  -  Mass  density  of  the  fluid. 

pc  -  Mass  density  per  unit  length  of  the  cable. 

o>  -  Angular  velocity  of  a  typical  cable  segment. 


CABLE  MODELLING 


Consider  a  finite  segment  model  of  a  cable  as  shown  in  Figure  2. 

The  segments  are  rigid  cylinders  connected  at  their  ends  with  spherical 
pins.  The  model  is  thus  a  chain  simulation  of  the  cable.  As  mentioned 
earlier,  recent  developments  in  the  analysis  of  multibody  system 
dynamics  has  led  to  computer  formulations  of  the  governing  equations 
of  motion  of  such  chain  systems.  These  formulations  have  the  option 
of  accomodating  models  with  segments  of  different  and  varying  lengths. 
The  shorter  segments  can  then  be  used  to  model  the  cable  in  those 
regions  where  the  cable  has  a  small  radius  of  curvature.  This  can  be 
accomplished  by  computing  the  angle  between  the  cylinder  axis  of  the 
segment  and  the  direction  of  the  resultant  force  vector  transmitted 
across  the  segment  joint.  If  this  angle  exceeds  an  arbitrarily  selected 
value,  the  segment  length  can  be  reduced.  (Recall  that  for  a  light 
flexible  cable,  the  resultant  force  at  a  cross  section  is  nearly 
tangent  to  the  cable.)  This  scheme  can  also  be  used  to  optimize  the 
number  of  segments  in  the  model. 


FLUID  FORCES  ON  THE  CYLINDRICAL  SEGMENTS 


Consider  a  typical  segment  of  the  cable  model  as  shown  in  Figure 
3,  Using  results  of  Webster  [35],  the  fluid  force  f  per  unit  length 
at  a  typical  point  P  may  be  expressed  as: 


f  =  a  aN  +  B!yN|vN  +  ci yT! vT 


(1) 


where  is  the  normal  component  of  the  fluid  velocity  relative  to 
the  cable  segment  at  P,  a^  is  the  normal  comoonent  of  the  fluid  accelera¬ 
tion  relative  to  the  cable  segment  at  P,  VT  is  the  tangential  component 
of  the  fluid  velocity  relative  to  the  cable  segment  at  P,  and  the  coef¬ 
ficients  A,  B,  and  C  are: 


A  «  CMp(ff/4)d2 

(2) 

B  *  CNp(d/2) 

(3) 

and 


C  -  CTPu(d/2)  (4) 

where  p  is  the  fluid  mass  density  and  d  is  the  diameter  of  the  cylindrical 
segment.  CM,  C^,  and  Cj  are  coefficients  dependent  upon  the  Reynolds 
number  of  the  fluid  flow  past  the  cable  segment.  These  are  usually 
determined  experimentally  and  the  reported  results  may  vary  slightly. 
Webster  I35J,  for  example,  records  them  as: 


( 0.0  for  R  <  0.1 
e  — 

0.45  +  5.93/(Re)0*33  for  0.1  <  R  <400. 

'  5 
1.27  for  400.  <  Rg  <  103 

0.3  for  R  >  105 

^  e  w 


(5) 


(6) 


and 


1 .88/(ReT)°* 74  for  0.1  <  RgT  <_  100.55 
0.062  for  RgT  >  100.55 


(7) 


where  the  Reynolds  numbers  RgN  and  RgT  are  defined  as 


Re  a  pd|VN|/u 


(8) 


and 


Rei  =  pd|VT[/y 

where  u  is  the  viscosity  of  the  fluid. 


(9) 


Of  these  coefficients  CN>  the  normal  drag  coefficient  is  undoubt¬ 
edly  the  most  significant  in  affecting  cable  motion.  It  is  also  the 
most  widely  studied  in  the  literature  (See  for  example.  References 
(47-54].)  The  variation  of  CN  with  RgN  is  often  depicted  graphically  as 
in  Figure  4.  A  linear  curve  fit  of  this  graph  leads  to  the  expressions: 


f 1  *45  +  8.55  (Rg)"0'40  for  1  <  Rg 

1.0  +  4(R J'1/2  for  30  <  R  <  IOC 
e  —  e 

2.25  -  0.45  1  og( R  )  for  100  <_  R 


CN  =  < 


0.90  for  1000  <  R  <  4000 
—  e 

1.05  +  0.54  log(Re)  for  4000  _<  R 

1.21  for  15,000  <  R  <  150,000 
—  e 

,  0.3  for  150,000  <  R 

^  -  e  a 


<  30 


<  1000 

<  15000 


(10) 
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Figure  4.  Normal  Drag  Coefficient  (References  [47-54]) 


As  mentioned  aoove,  these  expressions  are  of  slightly  different  form 
than  those  in  Equation  (6). 

In  Equation  (1),  the  first  term  is  due  to  the  acceleration  of  the 
fluid  relative  to  the  cable  segment.  It  is  sometimes  called  the  "added 
mass  force"  (See  References  [55-571).  The  second  term  is  the  "drag  force. 
The  third  term,  the  "tangential  drag  force,"  is  tangent  to  the  cable 
segment  and  it  is  generally  smaller  than  the  other  fluid  forces.  The 
effects  of  each  of  these  forces  is  examined  in  more  detail  in  the 
following  sections  of  the  report.  It  is  assumed  throughout  the  analysis 
that  the  coefficients  A,  B,  and  C  of  Equation  (1)  are  constant,  or  at 
least  "slowly  varying"  along  the  length  of  the  cable  segment.  Also, 
the  fluid  is  assumed  to  be  incompressible. 


ADDED  MASS  FORCES 


RELATIVE  ACCELERATION  OF  THE  FLUID  NORMAL  TO  THE  CABLE 


To  determine  the  force  system  exerted  on  a  cable  segment  due  to  the 
added  mass  forces,  it  is  necessary  to  develop  a  description  of  aN,  the 
normal  component  of  the  fluid  acceleration  relative  to  the  cable  segment 
in  terms  of  the  ambient  motion  of  the  fluid  and  in  terms  of  the  motion 
of  the  cable  segment.  To  do  this,  let  G  be  the  mass  center  of  a  typical 
cable  segment  as  shown  in  Figure  3.  Let  x  locate  a  typical  point  P 
relative  to  G  and  let  L  be  the  length  of  the  cable  segment.  Finally,  let 
n  be  a  unit  vector  parallel  to  the  axis  of  the  segment  as  shown.  Then 
x  is  positive  if  the  vector  GP  has  the  same  sense  as  n. 


Let  aw  represent  the  acceleration  of  the  fluid  in  an  inertial 
reference  frame  R.  (Usually  §w  will  be  zero.)  Then,  at  point  P  the 
acceleration  of  the  fluid  relative  to  the  cable  segment  may  be  written  as: 


aw/P  „ 


(ID 


p 

where  a  is  the  acceleration  of  P  in  R.  ap  may  be  expressed  in  terms  of 
a^,  the  acceleration  of  G  in  R,  as: 


a6  *  ap/G 


(12) 


where  aP/^  Is  the  acceleration  of  P  relative  to  G  in  R.  However,  since 

P/G 

both  P  and  G  are  fixed  on  the  cable  segment,  a  may  in  turn  be  written  as 


(13) 


p /Q 

a  ~  aX(xn)  +  ux[wx(xn)] 

where  a  and  <o  are  the  angular  velocity  and  the  angular  acceleration 
of  the  segment  in  R.  Since  a^,  in  Equation  (1),  is  the  component  of  aW^ 
perpendicular  to  n,  it  may  be  written  as: 


*N  S 


.W/P 


-  (a 


W/P 


n)n 


(14) 


By  substituting  Equation  (13)  into  (12),  (12)  into  (11),  and  finally 
(11)  into  (14),  aN  becomes 

a  =  aW/G  -  xaXn  -  xu>X(uXn)  -  (aW/G  •  n)n  +  x(«  •  n)2n  -  xw2n  (15) 
-N  •  •  ~  ~  -  -  - 

where  a^G  is  aW  -  aG,  the  acceleration  of  the  fluid  relative  to  the 
cable  at  6  in  R,  and  u)2  is  a)  •  a).  Let  §GN  be  the  acceleration  in  R  of 
this  fluid  relative  to  G  perpendicular  to  the  segment  axis.  Then,  aGN  is 


aW/G 


(16) 


Hence,  from  Equations  (15)  and  (16),  aN  may  be  written  as 


a„  3  a-.,  -  xaXn  -  xwX(u>Xn)  +  x(ui 
-N  -GN  -  -  -  - 


xa>2n 


ft 


(17) 


RESULTANT  ADDED  MASS  FORCE  AND  MOMENT 


The  added  mass  force  at  a  typical  point  P  of  the  cable  segment  is 
AaN.  Let  the  set  of  added  mass  forces  acting  on  the  entire  cable  segment 
be  represented  by  an  equivalent  force  system  consisting  of  a  single  force 
F^  passing  through  G  together  with  a  couple  with  torque  T^.  Then  FM  and 
Tm  are  given  by 


L/2 

Em  =  /L/i/Vx 


L/2 

T  =  /  xnxAavdx  (19) 

-M  -N 

By  assuming  that  A  is  independent  of  position  along  the  cable  segment 
and  by  substituting  from  Equation  (17)  into  Equations  (13)  and  (19), 
Fm  and  Tm  become,  upon  integration 


-M  *  AL?GN 


T^  *  -A(L  /1 2)  (a  -  (a  •  n)n  +  (u  •  n)nXu»] 


n 


DRAG  FORCES 


RELATIVE  VELOCITY  OF  THE  FLUID  NORMAL  AND  TANGENT  TO  THE  CABLE  SEGMENT 


To  determine  the  drag  forces  exerted  by  the  fluid  on  the  model 
segments,  it  is  necessary,  from  Equation  (1),  to  have  expressions  for  V^ 
and  VT,  the  normal  and  tangential  components  of  the  fluid  velocity 
relative  to  the  segment  at  a  typical  point  P.  These  velocities  may  be 
obtained  by  using  the  same  procedures  outlined  in  the  foregoing  section 
of  the  report.  Specifically,  let  V  represent  the  fluid  velocity  in  R. 
Then,  at  P  the  velocity  of  the  fluid  relative  to  the  segment  may 
be  expressed  as 

vw/p  *  vw  -  vp  (22) 

P  r 

where  V  is  the  velocity  of  P  in  R.  Vp  may  be  expressed  in  terms  of  V  , 

the  velocity  of  G  in  R  as 

vp  *  yG  +  vp/G  (23) 

P/G 

where  V  is  the  velocity  of  P  relative  to  G  in  R.  However,  since  P 
and  G  are  both  fixed  on  the  cable  segment,  VP//G  may  be  written  as 

VP/G  *  uXxn  (24) 

where,  as  before,  w  is  the  angular  velocity  of  the  cable  segment  in  R. 

W/P 

Since  V^  is  the  component  of  V  perpendicular  to  n,  it  may  be  written  as 


V  =  vw/p 
:n  : 


(V 


,W/P 


n)n 


(25) 


By  substituting  Equation  (24)  into  (23),  (23)  into  (22),  and  finally 
(22)  into  (25),  VN  becomes 


VN  =  vW/G  -  (VW/G  •  n)n  -  uXxn  (26) 

where  V^G  is  VW  -  VG,  the  velocity  of  the  fluid  relative  to  G  in  R. 
Let  V™  be  the  velocity  of  the  fluid  relative  to  G  in  R  perpendicular 
to  the  cable  segment  axis.  Then  VGN  is 


(27) 


Hence,  from  Equations  (26)  and  (27),  VN  finally  becomes 


-N  =  -GN  "  -Xx- 


(28) 


Similarly,  since  Vy  is  the  component  of  parallel  to  n,  it 
may  be  written  as: 


VT  =  (VW/P  •  n)n 


(29) 


Using  Equations  (22),  (23),  and  (24),  this  may  be  written  as 


vT  *  (vw/G  •  n)n  =  vGT 


(30) 


n 


where  VgT  is  the  velocity  of  the  fluid  relative  to  6  in  R  parallel  to 
the  cable  segment  axis. 

Finally,  in  Equation  (1),  it  is  necessary  to  have  expressions  for 
the  magnitudes  of  VN  and  V-p  From  Equation  (28)  the  magnitude  of  VN 
may  be  written  in  the  form 

|VN|  =  (a  +  bx  +  ex2)1*  (31) 

where  a,  b,  and  c  are 

3  =  «GN  -GN  =  (-GN^2 

b  =  -2VGN  •  uxn  (33) 

and 

2 

c  =  (uxn)  •  (uxn)  =  (uxn)  (34) 

Similarly,  from  Equation  (30),  the  magnitude  of  Vy  may  be  written  as 


NORMAL  CRAG  FORCE  AND  MOMENT 


The  normal  drag  force  at  a  typical  point  P  of  the  cable  segment  is 
B|VN|VN  Let  the  set  of  normal  drag  forces  acting  on  the  entire  cable 
segment  be  represented  by  an  equivalent  force  system  consisting  of  a 
single  force  passing  through  G  together  with  a  couple  with  torque 
T^j.  Then,  F^  and  T^  are  given  by: 

L/2 

fn  =  b  /  |vN|vN  dx  (36) 

-L/2 


L/2 

TN  =  B  /  |VN|xnXyN  dx  (37) 

-L/2 

Consider  first  F^.  From  Equations  (28)  and  (31),  F^  may  be 
written  as: 

L/2 

FN  3  B  /  (a  +  bx  +  cx2)1*  (VGN  -  u>Xxn)dx  (38) 

-L/2 


L/2  L/2 

-N  *  ®-GN  f  (a  +  bx  +  cx2)ljdx  -  Buxn  /  x(a  +  bx  +  cx2)igdx  (39) 
-L/2  -L/2 

The  integrals  in  Equation  (39)  may  be  evaluated  in  closed  form.  In 
[57]  it  is  found  that  by  letting  X  be  defined  as 


(40) 


r 


x 


a  +  bx  +  cx 


2 


the  first  integral  may  be  evaluated  as 


L/2 

/  X^dx  =  [ (cL  +  b)/4c]Xh|  +  [(cL  -  b)/4c] X4] 
-L/2  L/2  -L/2 


+  [(4ac 


b2)/8c3/2]log(- 


X*5! 


.  LcJ  +  b 

LIZ _ t-SL.) 

lcH  .  J.  1 
L/2  7~  2c5 


(41) 


Similarly,  the  second  integral  in  Equation  (39)  may  be  evaluated  as: 


L/2 


/  xX^dx  =  (1/3c)X3/2|  -  (l/3c)X3/2|  -  [b(cL  +  b)/8c2] X^] 
-L/2  L/2  -L/2  L/2 


2 1  vh\ 


+  [b(b  -  cD/ScHX'  , 

-L/2 

-  [b(4ac  -  b2)/16c5/2]log(— ^ \ - ^L_) 

vHi  Lc 3  ,  b 
a  il/2-  -y  2cv2 


(42) 


Finally,  for  the  purpose  of  evaluating  T^,  it  is  useful  to  also  record 
the  result: 


L/2 

/  x2XJsdx  =  ( 1  / 4c )  [(L/2)  -  (5b/6c)]X3/2|  +  (l/4c)  [(L/2)  +  (5b/6c)]  r/l\ 

-L/2  L/2 

+  [(5b2  -  4ac)/16c2]  { [ ( cL  +  b)/4c]  X^l  +  [(cL  -  b)/4c]XJ|( 

L/2  .  -L/2 


Lc*5  ,  b 


+  [(4ac  -  b2)/8c3/2]log(-r 


L/2*  - *  2P? 

•L/2  T  ^ 


■)) 


Hence,  from  Equations  (39),  (41),  and  (42),  F^  may  be  written  as 


L/2 


(43) 


Fm  =  B{ojxn(1 /3c)  (X3/2 1  -  X3/2i  )[(b/2c)u»xn 

-N  -  -  -L/2  L/2 


cL  +  b\ *  /CL  -  b ^ y ‘2 1 


+  V  ^  4c  _L/2 

2  X\/2  +  H~+  )y 

V2  ,'^\h<  Lc~‘»  ,  b 

A  i|_/2  ‘  2  2c^ 


,  , 4ac  -  b  , 

+  (  775  )^°9  v  • 

8c 


(44) 


Next,  consider  T^.  From  Equations  (37),  (28),  (31),  and  (40), 
Tn  may  be  written  as: 


L/2 

Tn  =  B  /  xX*5  nX(VGN  -  uiXxn)dx 
-L/2 


(45) 


By  expanding  the  vector  product  becomes: 


L/2 


L/2 


L.  =  BnxVr,,  /  xX^dx  -  8I«  -  (wn)nj  /  x2X^dx  (46) 


-GN 


■L/2 


-L/2 


By  using  Equations  (42)  and  (43),  TN  takes  the  form 


IN  =  BnxVGN[(l/3c)X3/2|  -  (1/3c)X3/2 


L/2 


-L/2 


b(cL  +  b)  x4j  +  b(-cL  +  b) 

8c2  L/2  8?  -L/2 

?  y  |  t.  Lc^  +  b 

.  b(4ac  -ft  l0?  ;X|l^  ,2,  &  )) 


16c 


xl 


Lc*  '  K 
L/2  '  2  2c' 


+  B[(w*n)n  -  u] { (1  /4c)  [(L/2)  -  (5b/6c)]X  /  1  j_/ 2 


(47) 


+  (l/4c) [(L/2)  +  (5b/4c)]X3/2|  +  [(5b2  -  4ac)/16c2] 


-L/2 


L/2 


Lc^  b. 


+  ^  jh  _+  log  (  t,,L/'  lcT--F-)]} 


cL-b  J*. 

4c  -L/2  '  8c3/2' 


'l/2  *  ~r  * 


X\/2  “ 


2  2c 
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TANGENTIAL  DRAG  FORCE  AND  MOMENT 


The  tangential  drag  force  at  a  typical  point  P  of  the  cable  segment 
is  given  by  the  third  term  in  Equation  (1)  as:  C|Vy|Vp  Hence,  let  the 
set  of  tangential  drag  forces  acting  on  the  entire  cable  segment  be 
replaced  by  an  equivalent  force  system  consisting  of  a  single  force  Fy 
passing  through  G  together  with  a  couple  with  torque  Jp  Then  Fy  and 
Ty  are  given  by 

L/2 

F  =  C  /  |VJV  dx  (48) 

-i  _L/2  -i  -i 

and 

L/2 

Tt  =  C  /  | VT | xnxVT  dx  (49) 

-L/2  *  ~l 

Consider  first  Fp  By  substituting  for  Vy  from  Equation  (30), 

Fy  becomes 

El  ■  CL|Vgt|Vgt  (50) 

Consider  next  Tp  From  Equation  (30)  it  is  seen  that  Vy  is 
parallel  to  n.  Hence,  the  vector  product  is  zero  and  therefore  Ty  becomes 

IT  ■  0  (51) 


WEIGHT  AND  BOUYANCY  FORCES 


Consider  finally  the  weight  and  bouyancy  forces  acting  on  the  cable 
segment.  The  weight  forces*  may  be  represented  simply  by  the  single 
vertical  (downward)  force  W  passing  through  the  mass  center  G  as  shown  in 
Figure  5.  If  pc  is  the  mass  density  per  unit  length  of  the  cable  segment, 
then  W  may  be  expressed  as 

W  =  -PcgLk  (52) 

where  k  is  a  vertical  unit  vector  and  g  is  the  gravity  constant. 

The  bouyancy  forces  are  not  quite  as  simple  as  the  weight  forces. 

The  bouyancy  forces  are  due  to  the  hydrodynamic  pressure  forces  exerted 
on  the  cable  segment.  These  forces  are  normal  to  the  surface  of  the 
segment  and  they  increase  in  magnitude  linearly  with  the  depth.  For  a 
submerged  cylinder,  the  bouyancy  forces  may  be  represented  by  the  single 
vertical  force  B  passing  through  the  geometric  center  G  of  the  cylinder 
as  shown  in  Figure  5.  However,  since  the  cable  segment  represents  a 
finite  segment  modelling  of  a  continuous  cable,  the  ends  1  and  2  of 
the  segment  are  not  exposed  to  the  fluid.  That  is,  they  are  "shielded" 
from  the  fluid  forces  by  the  adjacent  joining  cable  segments.  Therefore, 
the  hydrostatic  forces  on  the  cable  segment  are  all  normal  to  the  segment 
axis.  If  the  force  B  is  represented  by  components  and  perpendicular 

*These  are  the  only  non-fluid  forces  considered  in  the  analysis.  They 
are  included  in  the  analysis  because  of  their  close  association  with 
bouyancy  forces. 
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and  parallel  to  the  cable  segment,  then  the  resultant,  equivalent  hydro¬ 
static  force  on  the  cable  segment  may  be  represented  by  the  single 
force  Bn  passing  through  G.  (This  is  also  seen  by  noting  that  the 
hydrostatic  forces  on  the  ends  of  a  submerged  right  circular  cylinder 
are  parallel  to  the  axis  of  the  cylinder  and,  hence,  may  be  represented 
by  the  component  By.)  The  vertical  force  B  of  Figure  5.  may  be  expressed 
in  the  form  (See,  for  example.  Reference  [49].) 


B  *  -pgVk 


(53) 


where  p  is  the  mass  density  of  the  fluid  and  V  is  the  volume  of  the 
submerged  segment.  Hence,  since  B^  is  the  component  of  B  perpendicular 
to  the  segment  axis,  it  may  be  represented  as: 


?N  =  -  "  (-  *  n^n  =  Ox^xn) 


(54) 


or 


=  -pgnx(kxn) 


(55) 


on 


CONCLUSIONS 


To  summarize  the  foregoing  results,  the  force  system  exerted  on  a 
segment  of  a  finite  segment  cable  model  by  the  fluid  inertia,  drag,  and 
hydrostatic  forces,  together  with  the  gravitational  forces,  may  be 
represented  by  a  single  force  F  passing  through  the  mass  center  of  the 
cable  segment  together  with  a  couple  with  torque  T,.  F  and  T  are 
given  by  the  expressions 

-  =  -M  +  -N  +  -T  +  B~n  +  -  (56) 

and 


where  FM,  F^,  F^.,  B^,  W,  TM,  and  are  given  by  Equations  (20),  (44), 
(50),  (55),  (52),  (21),  and  (47)  respectively. 


Equations  (56)  and  (57)  are  in  a  format  which  may  be  directly  converted 
into  subroutine  algorithms  for  use  in  finite  segment  computer  codes 
(such  as  UCIN-SUPER  (46])  written  specifically  for  multibody  dynamics 
analyses.  That  is,  given  suitable  initial  conditions  for  the  kinematics 
of  the  cable  segments  together  with  the  segment  physical  and  geometrical 
parameters,  and  the  fluid  properties,  the  Reynolds  numbers  and  drag 
coefficients  can  be  determined.  This,  in  turn,  determines  F  and  T 
which  may  then  be  used  in  the  governing  equation  algorithm.  Upon 


numerical  integration,  the  kinematics  at  the  end  of  a  short  time 
interval  is  determined.  The  above  process  can  then  be  successively 
repeated  until  the  time  history  of  the  cable  motion  and  dynamics  is 
known. 

Reports  on  the  development,  validation,  and  application  of  such 
computer  algorithms  and  codes  are  currently  being  planned  and  prepared. 
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